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1.  INTRODUCTION 


This  paper  is  concerned  with  a  class  of  linear  predictor-corrector 
digital  filters  of  the  third  order  with  three  parameters.  This  class 
includes  the  four  digital  filters  that  have  been  or  are  being  used  at 
White  Sands  Missile  Range  (WSMR)  of  the  United  States  Army.  Two  of  them, 
namely  the  McCool's  QD  filter  [3]  and  Shepherd's  third  order  parabolic 
spline  filter  [A],  were  devised  at  WSMR  for  missile  trajectory  smoothing, 
reconstruction  and  differentiation.  All  of  them  have  applications  in 
real  time  radar  tracking  digital  servomechanisms.  The  theory  for  this 
paper  was  developed  largely  from  internal  notes  of  WSMR  by  Chui  [1]  and 
Shepherd  [A]  although  some  laborious  detail  in  [1]  and  [A]  is  not  given 
here.  The  class  of  linear  predictor-corrector  digital  filters  to  be 
discussed  here  can  be  defined  as  in  the  following.  Set 

w  -  (Wj,  w2,  w3]T  .  (1.1) 

Here,  and  throughout,  the  superscript  T  indicates  the  transpose  of  a 
matrix.  Hence,  w  in  (1.1)  is  a  three-dimensional  constant  column 
vector.  If  {x^},  1  <*  0,  1,  . ..,  denotes  an  input  signal,  h  a 

positive  constant,  and  the  sequence  of  three-dimensional  vectors 

*i  “  lyi*  yI’  yi,T»  (1.2) 

i  =  0,  1 .  denotes  the  output  response  corresponding  to  the  input 

signal  (x^J  and  subject  to  the  initial  condition  *  (y_i»  yIi»  y"j_)  * 

then  the  class  of  digital  filters  we  consider  is  defined  by 


If  che  response  due  Co  arbitrary  initial  value  is  to  damp  out,  it 

is  apparent  from  (1.7)  that 


lim  Anx  =  0  for  all  x  in  F^. 

n-x» 


0.8) 


The  relationship  between  this  observation  and  the  stability  of  the  filter 

(1.5) ,  (1.6)  is  developed  in  section  2.  More  precisely,  we  will  show  that 
if  all  the  eigenvalues  1^,  l2»  ^3  of  the  matrix  A  lie  in  the  open 

unit  disc  |z|  <  1,  then  the  digital  filter  (1.3)  -  (1.4),  or  equivalently 

(1.5)  -  (1.6),  is  stable.  Furthermore,  the  condition  | A ^ |  <  1,  for 

j  =  1,  2,  3,  is  satisfied  if  and  only  if  the  condition  (1.8)  is  satisfied. 

We  will  study  this  eigenvalue  problem  via  z-transforms.  This  technique 

will  be  extended  to  "uncouple"  the  system  (1.3)  -  (1.4)  solving  for  each 

of  the  output  response  sequences  (y^),  fypK  80(1  in  terms  of  the 

input  signal  (x^)  and  the  initial  conditions  y_^,  yV^,  and  y"^ 

respectively  in  recursive  forms.  These  three  recursive  formulas  will  allow 

the  studying  of  each  of  the  three  output* response  sequences  {y^},  ^yp^» 

and  (y")  individually  without  referring  to  the  other  two.  In  particular, 
P 

design  criteria  can  be  formulated. 

In  section  3,  we  will  consider  the  special  case  of  Shepherd's 
"parabolic  spline  predictor  corrector  filter",  which  is  obtained  from 

(1.5)  and  (1.6)  by  setting 


w  “ 


(1.9) 


where  6  is  considered  as  a  design  parameter. 
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In  section  4  we  consider  somewhat  more  briefly  McCool's  "QD"  filter 
{3],  for  which  (with  w  =  (a,  g/h,  Y/h2]T) 

fa.a.Yl  -P - 3 - - »  6  =  IT  »  V  =  Tfl  ’ 

l10H  +  331T  +  23M  -  6  "  M  -l 

Morrison's  "fading  memory  polynomial  filter  of  degree  2"  [5],  for  which 

[a,e,Y]  =  [1- 03,  |(l-0)2(l  +  6),  (1-0)3J, 

and  an  a  -  8  filter  (for  which  one  can  set  w^  =  0  and  y"^  “0  in 
the  beginning)  studied  by  Gonzales  in  1968  (cf.  [2]). 


2.  STABILITY  AND  UNCOUPLING  OF  THE  FILTERS 


We  will  use  the  technique  of  z-transforms  to  study  the  class  of 
digital  filters  defined  by  (1.3)  -  (1.4),  or  equivalently  (1.5)  -  (1.6). 

If  {b^}  is  a  bi-infinite  sequence  of  complex  numbers,  then  the  z-transform 
of  the  sequence  {b^},  j  *  . . . ,  -1,  0,  1,  . . . ,  is  the  formal  Laurent 
series 


B(z)  =  l  b  z~^  . 


Let  (x  ),  (y  },  {y'l  and  {y"}  be  defined  as  in  section  1.  We  will 
P  P  P  p 

set  xp  =  0  if  p  <  0,  and  yp  -  y£  -  yp  -  0  if  p  <  -1.  Hence,  the 
z-transforms  of  these  sequences  are  given  by 


X  =  X(z)  =  I  x  z  p 


respectively.  If  we  re-write  the  matrix  system  (1.5)  in  the  form  of  a 
system  of  three  simultaneous  difference  equations  and  take  the  z-trans forms 
of  each  of  these  three  equations,  we  obtain  the  following  system  of 
simultaneous  linear  algebraic  equations 


(z-l+w^  -  (h  -  hw^)Y2  -  (|h2-  |h2w1)Y3  *  w 


W2Y1  +  <Z“l  +  hw2)Y2  -  (h--|h  w2)Y3  -  w2zX 


w^^  +  hw3*2  +  w3^Y3  =  w3zX  • 


(2.1) 


In  matrix  representation,  (2.1)  can  be  written  as 


(A  -  2l3)(Y1,  Y2,  Y3K  =  -zXw 


(2.2) 


where  I3  is  the  3x3  identity  matrix.  Let  H^z),  H2(z),  and  H3(z) 

be  the  transfer  functions  of  this  digital  filter;  that  is. 


Yj (z)  -  H  <s)X(x) 


(2.3) 


for  j  *  1,  2,  3.  Then  H^z),  H2(z),  and  H3(z)  can  be  obtained  by 

solving  the  linear  system  (2.2)  using  Cramer's  rule.  Hence,  they  are 
rational  functions  in  z  *  with  the  same  denominator  det(  A  -  zl3).  This 


shows  that  if  det(A  -  zl.j)  jt  0  for  all  z  with 


<  1  or  Izl  >  1, 


or  equivalently  all  the  eigenvalues  of  A  lie  inside  the  unit  circle  |z|  =  1, 
then  the  filter  (1.3)  -  (1.4)  is  stable.  Let  X^,  X2,  X3  be  the  (not 

necessarily  distinct)  eigenvalues  of  A  and  let  x^,  x^  be  three 

*1 

corresponding  linearly  independent  eigenvectors.  Then  for  any  x  (  B  , 
x  -  +  o2x2  +  <*3X3  for  some  constants  a^,  a2,  and  a3>  Hence,  for 

any  positive  integer  n,  we  have 


5 


.n  ,  n  .  ,n  ,  ,  n 

A  x  -  +  o^  . 


This  shows  that  if  |Aj  |  <1  for  j  =  1,  2,  3,  then  Allx  £  as  n  -*•  °°. 

Conversely,  if  Anx  -*■  ()  for  all  x  €  1R  ,  we  can  pick  x  c  x (j  =  1,  2,  3) 

to  yield  X”x.  =  Anx.  -*•  0  as  n  -*•  °°.  Since  x.  ^  0,  we  mast  have 
7  J-J  —}  ~  “J 

|XX  |.  IX2 1 ,  |X3(  <  1.  That  is,  we  have  established  the  following 


THEOREM  2.1.  Let  X^,  X2,  X3  be  the  eigenvalues  of  the  matrix  A  giver: 
in  (1.6)  and  A  *  max(|X3|,  1X2),  |X3|).  Then  the  digital  filter  (1.3)  - 
(1.4)  is  stable  provided  A  <  1.  Furthermore  this  stability  condition 


A  <  1  holds  if  and  only  if 


lim  A  x  *  0 


for  all  x  €  1R  . 


We  next  determine  the  transfer  functions  H^(z),  H2(z),  and  H3(z) . 

To  do  this,  we  first  compute  det(zl3  -A).  It  is  given  by 

det(z!3~A)  =  z^  +  {(w^-3)  +  w2h  +^w^h2}z2 

■  r /  ,  \  .  .  .  1.2,  .  ,  to 


+  {(3-2w3)  -  w2h  +  yw3h  +  («]_  “  1)  •  (2.4) 


Hence,  by  Cramer's  rule,  we  obtain 


w  (w1~l)h 


'iw  -  <. ,(;!*-*)  '2  *+“2h-1 


i(w3  -  l)h2 


1  .2  . 

-jw2h  -  h 


,  1  1  2  . 

z  +  2  w3h  -  1 


(2.5) 


dTt^-A')  {V  +  feW3h  +  W2h  ‘  few3h  'W2h  +W1 
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Similarly,  we  obtain 


Vz)  =  detfll*-"A)  (2‘1)(V  +  hw3  “  W2>  * 


Y3(Z)  “  det(zIj-A)  (z-1)  w3  * 


Hence,  the  transfer  functions  H^(z),  H2(z),  and  H^(z)  as 

defined  in  (2.3)  can  be  written  as 


W1  +  W2^  ~  2w^)z  ^  +  ~  w2^*  +  w^)z 


-1  .  fl  .  L  2 


H, (z)  = 


1  1  +  ((w^-3) +w2h +-|w3h2}z  ^  +  {(3-2Wj)  -  w^h +^w3h^}z  2  +  (w-l)z  ^ 


and 


H2(z) 


w2  +  (hWj  -  2w2)z  1  +  (w2  -  hw^)z  2 


l+tfa^-S)  +  w2h  +^w^h2}z  1  +  {(3-2w1>  -  w2h +-^w^hZ}z  Z+(w1~l)z 


1  ,2,-2 


H3(z) 


,  -1  .  -2 
w^  -  2w^z  +  w^z 


1.  .2,-1 


1  +  ((wj-3)  +w2h  +“Wjht}z  A  +  {(3-2w1>  -  «2h +^w3h2}z  2  +  (Wj-l)z 


-3 


If  we  put  the  expressions  (2.8),  (2.9),  and  (2.10)  back  into  (2.3), 


multiply  the  denominator  of  H^(z)  to  *j(*)»  and  take  the  inverse  z-transforms 
of  each  of  the  expressions  for  j  ■  1,  2,  and  3,  we  obtain  the  following 


THEOREM  2.2.  The  digital  filter  given  by  (1.3)  -  (J.4)  can  be  written  as 
three  uncoupled  recursive  digital  filters: 

yp  ■  -{(Wj-3)  +w2h  +  2W3h2frp_i  “  {<3_2wi)  “w2h+Iw3h2^yp-2 

-  (w1-l)yp_3  +  Vp  ♦  (7'3h2+w2h-2w1)xp  l+  (|w;jh2-w2h+w1) 

yp  -  -{(Wj^-3) +w2h+^w3h2}y^_1  -  {(3-2w1>  -  w2h +|w3h2)y^_2 
~  («r1)yJ-3  +  W2Xp  +  (hw3-2w2)Xp-l  +  <W2'hw3)xp-2  ' 


(2.6) 

(2.7) 

(2.8) 

(2.9) 

<2.10) 


,(2.11) 

(2.12) 
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and 


y»  .  -Uwj-aj+Wjh+iwjh2},^  -  fo-awjj-Ujh+ivjh2},"  2 


-  (w  -l)yM  ,  +  w.x  -  2wnx  -  +  w^x  0  # 
1  p-3  3  p  3  p-1  3  p-2 


(2.13) 


where  p  =  0,  1,  2,  ...  with  initial  conditions  y  y'^,  and  y"  ,  and 
where  x  =  0  for  p-  <  0  and  y  =■=  y'  =  y"  =  0  for  p  <  -1. 

We  now  return  to  study  the  stability  of  the  filter  (1.3)  -  (1.4)  a 
little  closer  via  the  transfer  functions  H^z),  H2(z),  and  H3(z). 

From  (2.4),  putting  z  »  0,  we  see  that  det  A  =  1-w^  This  says  that 
A  =  0  is  an  eigenvalue  of  A  if  and  only  if  w  ■  1.  If  w3  =  0,  then 

the  transfer  function  H^z)  is  identically  zero.  On  the  other  hand, 

if  w 3  f  0,  then  the  value  z  ^  =  1  is  not  a  zero  of  H^fz),  while  it 

is  a  double  zero  of  H^z)  and  at  least  a  simple  zero  of  H2(z).  Hence, 

in  this  case,  all  the  eigenvalues  A^  A3  of  A  are  essential  in 

the  consideration  of  stability.  More  precisely,  we  have  the  following. 

COROLLARY  2.1.  Let  w3  j*  0.  Then  the  digital  filter  (1.3)  -  (1.4)  is 
stable  if  and  only  if  A  :=  max(|A1|,  |A2|,  |A2|)  <  1. 


3.  A  PARABOLIC  SPLINE  PREDICTOR-CORRECTOR  DIGITAL  FILTER 

|T 


Tliis  filter  is  the  special  case 


n  2(5 

«  “  \6,  T. 


26"] 

k2 


It  was  devised 


by  Shepherd  [4]  and  has  the  following  property:  Set 


vp(t)  ~yp  +  y’(t-ph)  +i[-^fyp-ify;+  <1-«>*;;+;TVi](t"ph>2’  (5 

1 

and  define  v(t)  on  [0,  “>)  by 

v(t)  ■  v  (t)  for  Ph  <  t  <  Ph  +  h,  p  =  0,  1,  2 .  /- 

p  —  V  ^ 

T 

Here,  the  input  signal  {x^}  and  the  output  response  =  [yp,  yp,  yp] 

satisfy  the  filter  relationship  (1.3)  -  (1.4)  with  w  =  (6,  -r~,  ~4,  J1. 

h  h 

Using  this  relationship,  it  is  easy  to  verify  that  the  function  v  defined 

by  (3.2)  is  indeed  a  parabolic  (or  third  order  C1)  spline  on  (0,  <»)  with 

knots  at  {0,  h,  2h,  ...}.  Since  v(ph)  =  y  and  v' (ph)  =  y',  the 

P  P 

spline  function  v  together  with  its  derivative  v'  interpolate  the 

T 

output  response  [y^,  yM  at  the  tine  t  *  ph,  p  =  0,  1,  ...  .  We 

like  to  think  of  the  input  signal  {x^}  as  a  noisy  measurement  of  (f(ph)) 

where  f  is  a  fairly  smooth  function  to  be  reconstructed,  and  the  output 

response  -  [yp,  ,  y " ]  as  an  approximation  to  [f(ph),  f ' (ph) ,  f"(ph)]  . 

The  approximation  improves  as  6  tends  to  1  as  can  be  seen  in  (3.3) 

below.  Hence,  the  spline  function  v  can  be  considered  as  an  approximation 

to  the  function  f  which  we  wish  to  reconstruct  from  the  noisy  data  (x  ). 

P 

This  particular  filter  is  a  one-parameter  digital  filter,  parametrized 

by  6  .  In  fact,  if  we  put  w  =  (5,  in  (1.3),  then  (1.3)  -  (1.4) 

It 

becomes 


Vl " 

(1  -  6)y] 

•  _ 

26 

yp+i 

h  yp 

m 

26 

"  h2^ 

rp  T  2' 


,»•  j.  x 

p  h  p+1  * 


25  25  ,  ,  ,,  ..  „  .  26 

-  77  »»  -Tyi  +  U-S)y 


"  h2  p+l  ’ 
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and  in  matrix  form,  we  have 


_  ,  tx  z6  25  ,T 

•2p+l  A6*p  +  Xp+1^’  h  ’  h2  1  »  wlth 


1-6 

(1  -  6)h 

|(1  -  6)h 

26 
~  h 

1-26 

(1-  6)h 

_26 

_  26 

1-6 

If  we  put  6-1,  then  y^+^  =  Xp+1  (3.3).  Hence,  if  6  is  close 

to  1,  the  spline  function  v  is  indeed  an  "approximation"  to  the  noisy 

input  data  {x^}.  However,  a  digital  filter  must  be  stable.  Therefore,  we 

will  discuss  how  close  can  6  approach  1  so  that  the  digital  filter 

2 

defined  by  (3.3)  remains  stable.  Since  w^  =  2 6/h  ^  0,  we  conclude  from 

Corollary  2.1  that  the  digital  filter  is  stable  if  and  only  if 

:=  max(  |  Aj,  | ,  |X2J,  |A^|)  <1  where  X^,,  X^,  and  X^  are  the  eigen¬ 
values  of  A^.  If  z  is  one  of  the  three  eigenvalues  A^,  X2,  and  X^, 
then  by  (2.4),  z  must  satisfy  the  equation 


z3  +  (46  -  3)z2  +  3(1  -  6)z  -  (1  -  6)  -  0. 


(3.5) 


This  equation  does  not  contain  the  time  increment  h.  Hence,  A^  is 
independent  of  h.  For  6*1,  the  solutions  of  (3.5)  are  0,  0,  -1, 
so  that  A^  =  1.  Hence,  the  filter  is  not  stable  if  6=1.  Let 

6  6n  "  1  “  n  ’  n  "  ±2’  *3’  •••  (3.6) 

so  that  6  -*•  1.  The  following  tables  indicate  the  stability  of  the  filter 
for  different  values  of  n.  Note  that  for  n  =  3,  the  roots  of  (3.5)  are 


i,  1,  -1,  so  that  A.  =  1  and  the  filter  Is  unstable.  Indeed,  Steelman  [5] 

O  63 

pointed  out  that  this  filter  Is  stable  If  and  only  if  2/3  <  6  <  1  although 
a  proof  Is  not  given  In  [5].  These  tables  verify  the  truth  of  the  statement. 
Note  that  In  Table  1A,  for  n  *  3,  6n  should  be  2/3  Instead  of  .667  and  the 
filter  should  be  unstable  as  mentioned  above. 


TABLE  1A.  STABILITY  VS.  MAXIMUM  MODULUS  OF  THE  EIGENVALUES  OF  A] 


_N 

DELTA  =  1  -  1/N 

MAXIMAL  MOD  OF  THE  ROOT 

STABILITY 

2 

0.500 

1.1228020 

NO 

3 

0.667 

0.9999995 

.NO 

4 

0.750 

0.9158811 

YES 

5 

0.800 

0.8543339 

YES 

6 

0.833 

0.8067909 

YES 

7 

0.857 

0.7685731 

YES 

8 

0.875 

0.7369311 

YES 

9 

0.889 

0.7101327 

YES 

10 

0.900 

0.6870247 

YES 

11 

0.909 

0.6668116 

YES 

12 

0.917 

0.6489261 

YES 

13 

0.923 

0.6329178 

YES 

14 

0.929 

0.6184893 

YES 

15 

0.933 

0.6054116 

YES 

16 

0.938 

0.5934409 

YES 

17 

0.941 

0.5824399 

YES 

18 

0.944 

0.5722803 

YES 

19 

0.947 

0.5627684 

YES 

20 

0.950 

0.5539731 

YES 

21 

0.952 

0.5458879 

YES 

22 

0.955 

0.5380396 

YES 

23 

0.957 

0.5307868 

YES 

24 

0.958 

0.5239464 

YES 

25 

0.960 

0.5174800 

YES 

26 

0.962 

0.5113528 

YES 

27 

0.963 

0.5055345 

YES 

28 

0.964 

0.5000002 

YES 

29 

0.966 

0.5835004 

YES 

30 

0.967 

0.6170683 

YES 

31 

0.968 

0.6419301 

YES 

32 

0.969 

0.6621327 

YES 

33 

0.970 

0.6793841 

YES 

34 

0.971 

0.6942258 

YES 

35 

0.971 

0.7074718 

YES 

36 

0.972 

0.7193651 

YES 

37 

0.973 

0.7301440 

YES 
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TABLE  1A  (Cont) 


_N 

DELTA  =  1  -  1/N 

MAXIMAL  MOD  OF  THE  ROOT 

STABILITY 

38 

0.974 

0.7399907 

YES 

39 

0.974 

0.7490339 

YES 

40 

0.975 

0.7573915 

YES 

41 

0.976 

0.7651453 

YES 

42 

0.976 

0.7723694 

YES 

43 

0.977 

0.7791147 

YES 

44 

0.977 

0. 7854404 

YES 

45 

0.978 

0.7913876 

YES 

46 

0.978 

0.7969885 

YES 

47 

0.979 

0.8022742 

YES 

48 

0.979 

0.8072767 

YES 

48 

0.980 

0.8120122 

YES 

50 

0.980 

0.8165183 

YES 

TABLE 

IB.  STABILITY  VS. 

MAXIMUM  MODULUS  OF  THE  EIGENVALUES  OF  A1 

_N 

DELTA  -  1  +  1/N 

MAXIMAL  MOD  OF  THE  ROOTS 

STABILITY 

2 

1.500 

3.4733100 

NO 

3 

1.333 

2.7423110 

NO 

4 

1.250 

2.3622870 

NO 

5 

1.200 

2.1263980 

NO 

6 

1.167 

1.9643850 

NO 

7 

1.143 

1  .-8455800 

NO 

8 

1.125 

1.7543650 

NO 

9 

1.111 

1.6819070 

NO 

10 

1.100 

1.6228290 

NO 

11 

1.091 

1.5736540 

NO 

12 

1.083 

1.5320190 

NO 

13 

1.077 

1.4962740 

NO 

14 

1.071 

1.4652290 

NO 

15 

1.067 

1.4379890 

NO 

16 

1.063 

1.4138780 

NO 

17 

1.059 

1.3923700 

NO 

18 

1.056 

1.3730690 

NO 

19 

1.053 

1.3556360 

NO 

20 

1.050 

1.3398040 

NO 

21 

1.048 

1.3253690 

NO 

22 

1.045 

1.3121390 

NO 

23 

1.043 

1.2999750 

NO 

24 

1.042 

1.2887430 

NO 

25 

1.040 

1.2783470 

NO 

26 

1.038 

1.2686840 

NO 

27 

1.037 

1.2596920 

NO 

28 

1.036 

1.2512910 

NO 

29 

1.034 

1.2434240 

NO 

30 

1.033 

1.2360490 

NO 
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TABLE  IB  (Cont) 


_N 

DELTA  -  1  + 

m 

MAXIMAL  MOD  OF  THE  ROOTS 

STABILITY 

31 

1.032 

1.2291180 

NO 

32 

1.031 

1.2225870 

NO 

33 

1.030 

1.2164250 

NO 

34 

1.029 

1.2105970 

NO 

35 

1.029 

1.2050840 

NO 

36 

1.028 

1.1998570 

NO 

37 

1.027 

1.1948880 

NO 

38 

1.026 

1.1901720 

NO 

39 

1.026 

1.1856750 

NO 

40 

1.025 

1.1813930 

NO 

41 

1.024 

1.1773080 

NO 

42 

1.024 

1.1734020 

NO 

43 

1.023 

1.1696650 

NO 

44 

1.023 

1.1660900 

NO 

45 

1.022 

1.1626620 

NO 

46 

1.022 

1.1593790 

NO 

47 

1.021 

1.1562260 

NO 

48 

1.021 

1.1531930 

NO 

49 

1.020 

1.1502780 

NO 

50 

1.020 

1.1474740 

NO 

NOTE: 

To  give  a  more 

accurate  picture ,  we  also  use  the  values 

6  =  N/SO. 

N  =  0, 

.  .  . ,  loo  as 

in  Table 

1C. 

TABLE 

1C.  Stability  VS. 

MAXIMUM  MODULUS  OF  THE  EIGENVALUES  OF  A] 

_N 

DELTA  =  N/50 

MAXIMAL  MOD  OF  THE  ROOT 

STABILITY 

0 

0.0000 

1.0040010 

NO 

1 

0.0200 

1.1561870 

NO 

2 

0.0400 

1.1873490 

NO 

3 

0.0600 

1.2053820 

NO 

4 

0.0800 

1.2170370 

.  NO 

5 

0.1000 

1.2247450 

NO 

6 

0.1200 

1.2296930 

NO 

7 

0.1400 

1.2325600 

NO 

8 

0.1600 

1.2337690 

NO 

9 

0.1800 

1.2336110 

NO 

10 

0.2000 

1.2325780 

NO 

11 

0.2200 

1.2299390 

NO 

12 

0.2400 

1.2266840 

NO 

13 

0.2600 

1.2226040 

NO 

14 

0.2800 

1.2177650 

NO 

15 

0.3000 

1.2122150 

NO 

16 

0.3200 

1.2060020 

NO 

TABLE  1C  (Cont) 


_N 

DELTA  =  M/50 

MAXIMAL  MOD  OF  THE  ROOT 

STABILITY 

17 

0.3400 

1.1991330 

NO 

18 

0.3600 

1.1916550 

NO 

19 

0.3800 

1.1835760 

NO 

20 

0.4000 

1.1749050 

NO 

21 

0.4200 

1.1656540 

NO 

22 

0.4400 

1.1558210 

NO 

23 

0.4600 

1.1454680 

NO 

24 

0.4800 

1.1344070 

NO 

25 

0.5000 

1.1227820 

NO 

26 

0.5200 

1.1105690 

NO 

27 

0.5400 

1.0977530 

NO 

28 

0.5600 

1.0842650 

NO 

29 

0.5800 

1.0700990 

NO 

30 

0.6000 

1.0552270 

NO 

31 

0.6200 

1.0396080 

NO 

32 

0.6400 

1.0232000 

NO 

33 

0.6600 

1.0059480 

NO 

34 

0.6800 

0.9877929 

YES 

35 

0.7000 

0.9686793 

YES 

36 

0.7200 

0.9484389 

YES 

37 

0.7400 

0.9270499 

YES 

38 

0.7600 

0.9043519 

YES 

39 

0.7800 

0.8801807 

YES 

40 

0.8000 

0.8543406 

YES 

41 

0.8200 

0.8265516 

YES 

42 

0.8400 

0.7964933 

YES 

43 

0.8600 

0.7638065 

YES 

44 

0.8800 

0.7275256 

YES 

45 

0.9000 

0.6870258 

YES 

46 

0.9200 

0.6407076 

YES 

47 

0.9400 

0.5859586 

YES 

48 

0.9600 

0.5176454 

YES 

49 

0.9800 

0.8165164 

YES 

50 

1.0000 

0.9999986 

NO 

51 

1.0200 

1.1474740 

NO 

52 

1.0400 

1.2783470 

•  NO 

53 

1.0600 

1.3992770 

NO 

54 

1.0800 

1.5134960 

NO 

55 

1.1000 

1.6228290 

NO 

56 

1.1200 

1.7284450 

NO 

57 

1.1400 

1.8311170 

NO 

58 

1.1600 

1.9314120 

NO 

59 

1.1800 

2.0297310 

NO 

60 

1.2000 

2.1263980 

NO 

61 

1.2200 

2.2216460 

NO 

62 

1.2400 

2.3156780 

NO 

63 

1.2600 

2.4086440 

NO 

64 

1.2800 

2.5006810 

NO 

65 

1.3000 

2.5918890 

NO 
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TABLE  1C  (Cont) 


_N 

DELTA  =  N/50 

MAXIMAL  MOD  OF  THE  ROOT 

STABILITY 

66 

1.3200 

2.6823650 

NO 

67 

1.3400 

2.7721790 

NO 

68 

1.3600 

2.8614030 

NO 

69 

1.3800 

2.9500870 

NO 

1.4000 

3.0382880 

NO 

71 

1.4200 

3.1260400 

NO 

72 

1.4400 

3.2133900 

NO 

73 

1.4600 

3.3003620 

NO 

74 

1.4800 

3.3869950 

NO 

75 

1.5000  ‘ 

3.4733050 

NO 

76 

1.5200 

3.5593280 

NO 

77 

1.5400 

3.6450760 

NO 

78 

1.5600 

3.7305680 

NO 

79 

1.5800 

3.8158280 

'NO 

80 

1.6000 

3.9008630 

NO 

81 

1.6200 

3.9856960 

NO 

82 

1.6400 

4.0703320 

NO 

83 

1.6600 

4.1547900 

NO 

84 

1.6800 

4.2390740 

NO 

85 

1.7000 

4.3232030 

NO 

86 

1.7200 

4.4071750 

NO 

87 

1.7400 

4.4910080 

NO 

88 

1.7600 

4.5747040 

NO 

89 

1.7800 

4.6582750 

NO 

90 

1.8000 

11.1541500 

NO 

91 

1.8200 

4.7706570 

NO 

92 

1.8400 

4.9082810 

NO 

93 

1.8600 

4.9914050 

NO 

94 

1.8800 

5.0744260 

NO 

95 

1.9000 

4.6608080 

NO 

96 

1.9200 

5.2401970 

NO 

97 

1.9400 

5.3229540 

NO 

98 

1.9600 

5.4056270 

NO 

99 

1.9800 

5.4882250 

NO 
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Hie  transfer  functions  Hfi(z)  :»  [H^z),  H2(z),  H3(z)]T  can  also  be 

of  or  ip 

obtained  by  substituting  w  -  (6,  -r-,  -f  ]  into  (2.8)  -  (2.10).  We  have 

n  a 


H  /_s  _  _ 6  +  6z _ 

”i'z'  -1  _2  -3  * 

1  +  (45-  3)z  +  3(1-  5)z  +  (5-  l)z 


li/i  -lx 
-r-(l  -  z  ) 

H  (z)  «  - -  - 

1  +  (46-3)z  A  +  3(1 -5)z  A  +  (5-l)z 


26,.  -1.2 
rjd  -  z  ) 

H  (z)  =  - ^ - - - —  , 

1  +  (46-  3) z  A  +  3(1  -  6)z  *  +  (6-  l)z  J 


Also,  by  applying  Theorem  2.2,  we  can  uncouple  the  system  (3.3)  to  yield: 

yp  m  ‘<46“3>yp_i  ■  3(1^5>yp-2  +  (1~  <S)yp_3  +  6xp  +  6*p_i  • 

,•  -  -(46  -  3),^  -  3(1  -  ♦  (1  -  «y'.3  ♦  f  *„  -  f  x 


26  46  .  26 


(3.7) 


(3.8) 


(3.9) 


(3.10) 


(3.11) 


yP  ”  -(46"3>y;-l  '  ^-^-2  +  (3-12) 


for  p  *  0,  1,  2,  ...  with  initial  conditions  x_^  >■  x  2  “  0  and 

T 

y_2  “  y_3  ”  o  and  _y_1  -  [y_lf  y^,  y^J  preassigned.  Equations  (3.10) 

(3.12)  enable  the  user  to  plot  the  output  response  *  (yp,  y^,  y ” J  in 

three  different  graphs,  where  the  graph  of  {yp}  shows  how  the  filter 

smooths  the  input  signal  (xp),  p  -  0,  1,  ...  .  It  is  also  interesting 

to  study  how  close  {y  }  is  to  (x  }  when  6  is  close  to  but  different 

P  P 

from  1.  This  can  be  done  using  the  transfer  function  H^(z)  given  in 

2 

(3.7)  as  in  the  following.  Since  the  filter  is  stable  for  <  6  <  1,  we 
always  pick  such  values  of  6.  For  aiaplicity,  we  set 


fit*,/'''  ' 
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e  *  1  -  6  >  0 


so  Chat  we  have 


1  -  H1(z) 


=  c 


-M 

(1-2  ) 


1  +  (l-4e)z  *  +  3ez  2  -  ez  ^ 


(3.13) 


By  applying  Parseval's  identity,  we  have 


,!i  w2  ■  &  C ,x<e^  ■  vri“)i2d" 

-  TH  f  I1  -  |2  |x(u>)  |2  da. 

'-IT  ^ 


where  we  have  used  the  common  notation  X(o»)  ■  X(e  lW)  and  x_^  *  0. 
Hence,  from  (3.13),  we  have 


9  •  ,  -2  r*  .  iu>  3 

I  y  I2  +  Y  |x  _y  I2  »  £_  - 1-1.--?  ) - 

-1‘  n  P  2ff  J  i  .  i2ta>  i3ta> 

p=0  r  r  J-n  l  +  (l-4e)e  +  3ee  -  ee 


(3.14) 


|X(o.)  |  do. 


«!(’ _ 

■ -If  (1  +  P  COSO)  + 


(1  -  costa))' 


2  2 
3c  cos  2b)  -  £  cos  3b>)  +  (p£sinb>  +  3esin2ta>  -  r.sin3u>) 


|X(w) |2  du,  (3.15) 


where  Pc  :■  1  -  4e  . 

The  last  expression  allows  us  to  design  Shepherd's  parabolic 'spline 

predictor-corrector  digital  filter  in  an  optimal  way.  Depending  on  the 

input  spectrum  X(b)),  one  can  pick  e  -  1-6  (numerically)  such  that  0  <  e  <  1/3 

and  such  that  the  last  expression  above  is  minimized.  In  particular,  if  we 

have  faith  on  the  data  (x  },  then  we  can  pick  e  >  0  very  close  to  zero. 

P 

This  gives 
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|y-il  +  I  K-»„l‘ 


p=0 


4e 

Tf 


r- 

J-Ti  1  + 


(1  -  COS  U))' 


2p  cos  a)  +  p 
e  c 


[ X (ai)  | 2  doi  (3.16) 


where  Pr  =  1  -  4e  as  above.  This  integral  can  be  used  to  design  the 
filter  efficiently. 


4.  THE  a-g-y  TRACKERS 


The  a-g-y  tracking  eq.uations  discussed  by  Steelman  [5]  can  be  obtained 

2  T 

from  (1.3)  -  (1.4)  by  putting  w  ■  [a,  g/h,  Y/h  ]  .  This  shows  that  each 
a-g-v  tracker  is  a  third  order  predictor-corrector  filter,  and  conversely 
(e.g.  the  filter  of  section  3  is  an  a-g-Y  tracker  with  g  *  Y  »  26).  From 
(2.8)  -  (2.10),  the  transfer  functions  of  this  filter  can  be  found  to  be 


Hx(z) 


a  +  (-2a+  g+-|y)z-1  +  (u- g+-ij'()z“‘ 

1  1  I2  Tj  (4. 1 ) 

1  +  (a+g  +  |y-3)z  x  +  (-2a- g+|Y+3)z  z  +  (a-l)z  5 


and 


H2(z) 


H3(z) 


1  +  (a+  g  +-|y -  3) z  1  +  (-2a-  g  +  |y  +  3)z~2  +  (a-  l)z 

X. 

2  2  2 
h  a  h 


~  (4.2) 


1  +  (a+  g  +  iy-  3)z  1  +  (-2a-  g+^y+  3)z  2  +  (a-  1); 


-3 


(4.3) 


These  equations  were  also  obtained  in  (5]  in  a  different  form.  However, 
using  the  formulation  in  (4.1)  -  (4.3),  we  can  immediately  uncouple  the 
filter  as  in  section  2,  yielding 


yp  -  -<ct+B  +  2V'3)yp-1  +  (Za+B-yY-  3)yp  2  -  (a- l)yp_3 
+  axp  -  (2a  -  g  -  ^Y)xp_1  +  (a  -  g  +  |y)xp_2  , 


(4.4) 


•  *$■'*#<* 
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(4.5) 


y*  =  _(a  +  0  +  iy_3)y'  +  (2a  +  0 - -  3)y'_ 2 


V  (2(3  Y)  ,  (P  Y) 

h  p  'h  hJ  p-1  Mi  hy  p-2 


b 


y"  =  -(a+3+|y- 3)y"  +  (2a+b-±Y-3)y"  2  - 


+  X  x  _  2  —  x  +  X  x 

h2  P  h2  P'1  h2  P-2  • 


(a  -  l)y"  , 
p-3 


(4.6) 


where  p  =  0,  1,  2,  ...  ,  y  =  y'  =  y"  =  0  if  p  <  -1  and  x ■  =  x  . ,  =  0. 

*  P  P  P  -1-2 

The  case  where  ct=*  1-  03,  B  »  |(1  -  0)2(]  +  0) ,  and  Y  =  (1  -  0)  3 
is  called  a  "fading  memory  polynomial  filter  of  degree  2 "  by  Morrison 
(cf.  [5]).  The  uncoupled  recursive  filters  (A. 4)  -  (4.6)  can  be  simplified 
to  be 

t 

y  =  30y  .  -  302y  ,  +  03y  ,  +  (1  -  03)x 

p  yp-l  yp-2  yn-3  p 


30(1  -  02)x  ,  +  302(1  -  0)x  „  , 

p-1  p-2 


y;  -  3%-i  -  *  ^-3 +  !«- 6,2(1 +<n 

-  £<1- 0)2(l  +  20)x  .  +^-(1-0)2(1  +  5O)x 
n  p-1  2  ii  p-2 


(4.7) 


(4.8) 


y;  -  30y’;-i  -  ™2y'i- 2  +  03yp-3  +  ^~(xp  -  2Vi +  xp-2k  l*-*) 


for  p  =  0,  1,  ...  with  y  ■  y'  -  y"  ■  0  for  p  <  -1  and  x  .  *  x  _  =  0. 

p  p  p 


h^j**'*  JMM*  i- 
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Note  that  the  feed-back  coefficients  are-  particularly  simple.  The  stability 
of  this  filter  is  also  particularly  easy  to  check.  In  fact,  it  is  stable  if 
and  only  if  | @ |  <1. 


As  another  special  case  of  the  a-$~Y  filter,  let  us  set  y  =  0  . 

This  is  the  so-called  a-3  tracker,  studied  in  WSMR  by  Gonzales  in  1968 
(cf.  [2]).  In  this  case,  the  factor  1  -  z  1  can  be  cancelled  in  each  of 
the  expressions  in  (4.1),  (4.2),  and  (4.3);  and  the  transfer  functions 
of  the  ot-3  tracker  simply  become 


and 


H^z) 


U2(z) 


_ a  -  (a-g)z  _ 

1  +  (a +  3  -  2)z-1  -  (a  -  1 )  z 


1  +  (a  +  3  -  2)z  1  -  (a-i)z 


(4.10) 


(4.11) 


H3(z)  =  0. 


(4.12) 


Since  the  denominator  is  a  quadratic  polynomial,  we  know  immediately 
that  this  a-3  filter  is  stable  if  and  only  if  A  „  <  1  where 

7  a, 3 


2 

lienee,  it  follows  easily  that  for  3  >  [  (o  +  3)/2]  ,  we  have  A  .  <  1 

O  jP 

if  and  only  if  a  >  0.  That  is,  a  stability  condition  for  the  a-B  tracker 
is : 
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(4.14) 


a  >  0,  (a  +  P)2  _<  4P  . 

Another  sufficient  condition  for  stability  of  the  a-P  filters  is 

P  >  0,  a  +  P  <  2.  (4.15) 

Of  course,  the  other  sufficient  condition  for  stability 

a  >  0,  P  >  0,  2a  +  P  <  4,  (4.16) 

which  also  follows  easily,  was  obtained  by  Gonzales  (cf.  [2,5]).  Using  an 
inverse  z-transform  as  in  section  2,  the  a-P  tracker  equations  can  be 
uncoupled  into  the  form: 

yp  =  -(a  +  P  -  2)yp  l  +  <°-l>yp_2  +  axp  ~  <“-0>xp-i*  (4.17) 

yp  ■  -<«  +6  -  2)yp_1  +  <«  -  i)yp_2  +  £<xp  -  xp_i>»  (4-18) 

and  y"  =  0,  p  =  0,  1,  ....  where  y  =  y'  =  0  if  p  <  -1  and  x  =0 
p  P  P  P 

if  p  <  0.  The  condition  (4.15)  is  very  useful.  It  says  that  if  the  first 
feed-back  coefficient  is  nonnegative  and  P  >  0,  then  the  a~P  filter 
(4.17)  -  (4.19)  is  always  stable. 

Finally,  let  us  discuss  McCool's  QD  filter  Icf.  3],  namely: 


60M2 


10M3  +  33M2  +  23M  -  6 


P  =  ~  and  y 


2a 


H2  ’ 


(4.19) 

(4.20) 
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where  M  is  a  natural  number.  The  transfer  functions  of  this  filter  can  :»■ 
obtained  by  substituting  (4.20)  into  (4.1)  -  (4.3)  and  the  uncoupled  recur¬ 
sive  input-output  relationship  can  be  obtained  from  (4.4)  -  (4.6).  We  nv>u 
study  the  stability  of  the  QD  filter.  Set 

f  ( z )  =  Z3  +  [(M+L)2,..,  -  3)z2  +  ((-:>>!’  -  2M+ 1  ),'>  +  3  |z  +  fM2,  .  1) 

with  P..  =  60/ (10M3  +  33M2  +  23M  -  6).  If  A  A„ ,  and  A,  are  the 
ft  1  z  3 

roots  of  f(z)  =  0  and  A^  =  max(|A^|,  |A^|,  lA.J),  then  by  Corollary 
2.1,  since  w^  =  Y/h~  =  2P^  /  0,  we  note  that  the  QD  digital  filter  is 
stable  if  and  only  if  A  <  1.  The  following  table  indicates  the 
stability  of  this  filter  for  M  =  0,  ±1,  ±2,  ...»  120  . 
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r 

TABLE  2. 

1  ■  . . . 

STABILITY  FOR  McCOOL'S  QD  FILTER 

M 

ro  =  pm 

MAXIMAL  MOD  OF  THE  ROOT  STABILITY  ] 

-20 

-0.0009 

1.2024780 

NO 

-19 

-0.0011 

1.2146490 

NO 

-18 

-0.0012 

1.2282810 

NO 

-17 

-0.0015 

1.2437600 

NO 

-16 

-0.0018 

1.2615480 

NO 

»  ; 

-15 

-0.0022 

1.2820790 

NO 

-14 

-0.0028 

1.3061450 

NO 

-13 

-0.0036 

1.3346260 

NO 

-12 

-0.0047 

1 . 3689580 

NO 

. 

-11 

-0.0063 

1.4111290 

NO 

-10 

-0.0087 

1.4641570 

NO 

-9 

-0.0124 

1 . 5328290 

NO 

* 

-8 

-0.0188 

1.6252010 

NO 

* 

-7 

-0.0303 

1.7560380 

NO 

-6 

-0.0538 

1.9554370 

NO 

-5 

-0.1099 

2.2953960 

NO 

i 

-4 

-0.2857 

3.0000000 

NO 

-3 

-2* 

-1.2500 

5.2570590 

NO 

-1 

-10.0000 

4.8910190 

NO 

0 

-10.0000 

13.5226500 

NO 

1 

1.0000 

0.9999990 

YES  , 

2 

0.2381 

0.5061458 

YES 

3 

0.0952 

0.6546553 

YES 

< 

4 

0.0478 

0.7329382 

YES 

5 

0.0275 

0.78191 32 

YES 

i 

6 

0.0172 

0.8155725 

YES 

k 

7 

0.0115 

0.8401833 

YES 

. 

8 

0.0081 

0.8589627 

YES 

I 

9 

0.0059 

0.8737922 

YES 

! 

10 

0.0044 

0.8857723 

YES 

i 

11 

0.0034 

0.8956649 

YES 

12 

0.0027 

0.9040021 

YES 

► 

13 

0.0022 

0.9110962 

YES 

14 

0.0018 

0.9172200 

YES 

■ 

15 

0.0014 

0.9225416 

YES 

; 

16 

0.0012 

0.9272274 

YES 

i  ; 

17 

0.0010 

0.9313622 

YES 

18 

0.0009 

0.9350671 

YES 

19 

0.0007 

0.9383648 

YES  ! 

20 

0.0006 

0.9413764 

YES 

i 

i 

★ 

Coefficients  of  the  polynomials  are  undefined;  division  by 

23 

zero. 
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